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en ! Abstract: We investigate cosmological consequences of nonlinear sigma model coupled 

> . 
^ . with a cosmological fluid which satisfies the continuity equation. The target space action is 



in 



of the de Sitter type and is composed of four scalar fields. The potential which is a function of 



^ ■ only one of the scalar fields is also introduced. We perform a general analysis of the ensuing 

cosmological equations and give various critical points and their properties. Then, we show 

CN ' that the model exhibits an exact cosmological solution which yields a transition from matter 

domination into dark energy epoch and compare it with the A-CDM behavior. Especially, we 

-^ , calculate the age of the Universe and show that it is consistent with the observational value if 

^ . theequationofthestatew/ofthecosmologicalfluidis within the range of 0.13 < uj < 0.22. 

c^ ■ 
" " " Some implication of this result is also discussed. 
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1. Introduction 

The recent astronomical measurements revealed that the current Universe is accelerating [1,2]. It is 
believed that the acceleration is caused by an unknown energy, i.e., dark energy, and grasping the identity 
of dark energy is one of most fundamental problems in the modem cosmology [3]. Many theoretical 
models for dark energy have been proposed ever since. Among them, the standard approach is to 
introduce a cosmological constant [4-9]. In spite of its simplicity and theoretical diversities, it confronts 
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the extreme fine tuning problem. An attractive alternative is to consider that the acceleration is driven 
by a scalar field. The well known models include quintessence [10,11] phantom [12], k-essence [13-15] 
with the non-canonical kinetic term for the scalar field, and quintom [16] models. In these models various 
aspects of dark energy can be well described in terms of dynamics of the scalar field and especially, the 
smallness of the cosmological constant is attributed to the decaying scalar energy density [17]. 

On the other hand, dark matter plays a central role in the early Universe in the process of structure 
formation [3]. Most of the dark matter is the cold dark matter (CDM) which is non-relativistic and 
non-baryonic. In particular, the CDM model with the cosmological constant A is established as 
a standard cosmological model, A-CDM [3]. This model is in good agreements with the CMBR 
data [18-20] as well as SN la data [1,2]. The special feature of A-CDM model is that it assumes the 
Universe with a zero spatial curvature. 

Recently, the dark energy model based on a nonlinear sigma model [21] with or without a 
cosmological constant was investigated. In this model the target space is noncompact four-dimensional 
de Sitter manifold and four scalar fields are introduced to account for this. It was found that an 
exponentially accelerating solution is possible even without the cosmological constant and that the model 
could describe dark energy interacting with stiff matter even without any matter present. The possibility 
of other kind of dark matter should be also addressed and it motivates to extend to the more general 
situation where matter is introduced separately. Therefore, in this paper, we consider more general 
case of the de Sitter nonlinear sigma model where the matter with the equation of state cuf and also an 
exponential potential term for the nonlinear sigma model are added to the action, and we investigate the 
cosmological consequences. 

Recall that in spite of the dynamical explanation of the smallness of the current dark energy density 
contrary to the cosmological constant, the approaches of [10-15] have their own shortcomings. For 
quintessence model, it has the nice properties of tracking solution and scaling behavior, but it is 
somewhat difficult to match the equation of state to be close to —1, and has its own fine-tuning 
problem [22]. The phantom model gives favorable result with the equation of state, but it has the issue 
of quantum instability, even though it is classically stable. The k-essence model has the ghost problem 
coming from the higher derivative nature of the theory. We do not attempt to address these difficult issues 
with the introduction of the four scalar fields, but one advantage is the existence of the exact solution, 
which makes comparisons with A-CDM more direct. ( See Equations (11) and (44).) 

Let us mention other salient features of our investigation. The field content consists of a phantom 
and the triplet fields which are canonical scalar fields. Unlike the ordinary phantom model where Big 
Rip singularity is known to exist, we have an explicit solution which does not show such a behavior. 
As far as the exact cosmic solution is concerned, the potential is of a negative exponential form of the 
scalar field 0. This solution prolongs from the matter-dominated epoch to the dark energy epoch. Recall 
that in the pure phantom model, the negative kinetic energy prohibits the potential to assume a negative 
value. In our approach, the triplet of scalar fields provides enough compensating positive energy density 
such that the weak energy condition is not violated. This exact solution can be exploited to extract 
some numerical values that can be compared with current observations. In particular, we calculate the 
age of the Universe and compare with the observations. We find that with a fine-tuning of a couple of 
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parameters, 0.13 < cuf < 0.22 is consistent with the observation. This matter corresponds to an exotic 
cosmological fluid. 

The paper is organized as follows. In Section II, we present a basic analysis of A— CDM model 
paying attention to the exact cosmological solution and its stability analysis. In Section HI, we consider 
the Einstein gravity coupled with the de Sitter nonlinear sigma model including an exponential potential 
and cosmological fluid with the equation of state cuf, and discuss the cosmological evolution equations. 
In Section IV, the stability analysis is performed and various critical points are identified. In Section V, 
we present the exact cosmological solution for a negative exponential potential and calculate the current 
age of the Universe. We compare it with the observational data and A— CDM model. Section VI includes 
conclusion and discussion. 

2. ACDM Model 

Let us first consider an action in which the Einstein gravity has a cosmological constant (c.c.) term 
with a matter term: 



■"matter 



(1) 



where k^ = SnG. Introducing the standard space-time metric via 

ds^ = -dt^ + a^{t)dxidx' (2) 

leads to the following equations: 

2 ^ ^ A 

3 3 ^ ^ 

H = - — {l + UJ.m)Pm (4) 

pm + 3{l+UJrn)Hpra = (5) 

In these expressions, Um is the equation of state parameter for the matter term, which satisfies 
Pm = ^mpm whcn T^ = {- pm, Pm, Pm, Pm) arc assumcd. 

In order to check the stability, we introduce the following dimensionless quantities (for pm, A > 0) 

K^/p^ y/X 

X = —7= — , y = —= — (6) 

From the above quantities, one obtains x' and y' as (N = In a) 

x' = ^ = l{l + um)x{x'-l) (7) 

The critical points, i.e., the solutions corresponding to x' = 0, y' = and the eigenvalues for the critical 
points [23] are given by 

3 

Xc = 0, yc = I {c.c. dominant), pi =0, P2 = — (1 + ^m) (9) 

3 

Xc = 1, yc = {matter dominant) , pi = 3{l + Um), P2 = 7:{^ + '^m) (10) 
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From these eigenvalues we see that the matter dominant phase is unstable for tUm. > — 1 while c.c. 
dominant phase has a decaying mode, i.e., 112 < being stable. (The existence of a zero eigenvalue 
(/ii = 0) in Equation (9) is originated from the fact that two variables x and y are connected by the 
relation x^ + y^ = 1. Therefore in this case one can reduce to one— dimensional space [24].) It is 
well-known that there exists the non-perturbative solution of Equations (3)-(5) which connects these 
two critical points as 

a(t) = a.(sinh[At])2/3a+--), H{t) = J ^coth[At], p^ = ^{smh[At])-^ (11) 

V o K 



where A = {1 + Wm) v 3A/2. This solution describes a smooth transition from matter-dominated power 
law expansion at early times into cosmological constant-dominated exponential acceleration at late times. 
To check the stability of the above solutions Equation (11), we consider the variation of 
Equations (3)-(5), which yields 



K^ 



2H5H = —5pm (12) 

6H = -y(l + w„)5p„ (13) 

Spm + 3(1 + UJm)p,nSH + 3(1 + UJm)H6pm = (14) 

Then we find the solution for the above Equations (12)- (14) as 

SH = BCtanh[At]e-f^'yK\ 6p„, = Ce'^'-'^ /k,^ (15) 

where B = {1+ u^)/AA, C is an arbitrary constant and f{t) = ~2At + ln[-l + e^^*] + 2 ln[tanh[At]]. 
Note that f{t) approaches 2 At as t — ?> oo and both 5H and 5pm decay, which implies that the 
solution (1 1) is stable. 

Introducing the dimensionless density parameters Vtm = K^pm/^H"^ and f^A = A/3iJ^, Equation (3) 
becomes 

fi„, + ^A = 1 (16) 

Particularly for the solution (11), Vt^ and f^A are given by 

^"^^ (cOsh[At])2' ^^ = ^- (cosh[At])2 ^^^^ 

By using the experimental data, i.e., px = ^ ~ (lO^^'^GevY ~ 10~^erg/c'm^, Hq ~ 2.28 x 10~^^s^^ 
and the solution Hq = -y3Acoth[ylto]/3, we can evaluate the value Ato = 1.27. From these values we 
find VIa and to for the dust- like matter (cUm = 0) as 

^A ~ 0.73, to ~ 4.34 X lO^^s (18) 

which are in agreement with the observational data [25]. 
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3. de Sitter Nonlinear Sigma Model with Potential 

In this section, we extend the analysis performed in the case of A— CDM to the de Sitter nonlinear 
sigma model [21] with a potential term. The cosmological fluid with the equation of state cuj is also 
added. The starting action is 

d'xy^ [ —R - ^G„^($)a^$"a,$^ - eV{(t>) + Cfiu,d ] (19) 

where $" = (</>, a*) {i = 1,2, 3), A is a dimensionless coupling constant and e = ±1. Here G^/^ is the 
metric of the de Sitter target space, 

G„^ = (-1, + e2<^, + e'^*^, + e'^*^ ) (20) 

where ^ is an arbitrary positive constant and V{(f)) is the potential given by 

¥{(!)) = Voexp{-K^(l)) (21) 

with an arbitrary constant 7 and Vq > 0. Among diverse possibilities, we have chosen exponential 
potential. A couple of reasons could be cited. The first one is that this potential is the prototype which 
gives rise to accelerating Universe and has interesting properties like scaling solution and attractor in the 
quintessence [26] or phantom model [23]. The second one is that it can yield non-perturbative solution 
of the type discussed in the previous section in our case for the negative potential with e = — 1. We 
regard the second and third terms of Equation ( 19) as representing the dark energy sector. For the matter 
part, we assume a cosmological fluid of the perfect fluid form, T^ = (— p/, Pf, Pf, pj), which satisfies 
the continuity equation, W ,_,T^^ = 0. 
We first note that the following ansatz 

a' = x' (22) 

solves the a* field equations. The above ansatz first appeared in higher dimensional gravity theory 
in association with spontaneous compactification of the extra dimensions [27-30]. It was revived in 
four dimensions recently in describing the accelerating Universe with the de Sitter nonlinear sigma 
model [21]. Note that it does not break the isotropy and homogeneity of the universe as long as we do 
not introduce the potential for the a fields. With this ansatz, the standard space-time metric Equation (2), 
and (f) = (l){t), the evolution equations, are given by 



H' 



H 



2^2 
3A2 

K 

~A2 



2 



2"^ 2K^a^ 2 



2 



l^' + Trie-^' + ^V^(0) 1 + ^Pf (23) 



3 



2 r 1 -, ,.2 






y(l + u;/)p/ (24) 



2kC0 fK'yX'^ 

= ^ + 3H<p-3^^^ + -i—Voe-^^'^ (25) 



and the continuity equation implies p/ = p(0)[a(t)/a(0)]^'^*^^+'^^'', where w/ is a barotropic equation of 
state with < w/ < 1. We mention a couple of properties of the above evolution equations. The first 
one is that in spite of the being a phantom, the presence of the second term in Equation (24) that comes 
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from the spatial variations of cr* fields could prevent the Big Rip singularity from being developed; it is 
not guaranteed that H will stay always positive at late times when pf is ignored. The other is that even in 
the e = — 1 case, the weak energy condition could not be violated (in the ordinary phantom model with 
a negative potential, the weak energy condition is always violated because of p<^ = —0^/2 — V{(j)) < 0). 
This again is due to the the second terms in Equations (23) and (24) coming from the spatial variations. 
In fact, we will show that for e = — 1 there exists an exact cosmological solution of the type discussed 
in A— CDM case, which interpolates between matter-dominated and dark energy-dominated epochs. 

4. Stability 

To perform the stability analysis and figure out the energy dominance of the kinetic (x), spatial (y) 
and potential (z) parts in Equation (23), we first introduce the following dimensionless quantities, 

X = —p^ , y = — -TT— , z = —= — (26) 

Then the constraint Equation (23) is given by 

-x' + y' + ez' + ^ = l ill) 

With A^ = In a, we obtain 



^ = -d^=2^-^^^-'^f^''' ~ ^^ + 3u;;)y^+3(-l+a;;) - H^+ujf)A ' ^z' +iVQy' (28) 
^ = I [-3(1 - ujf)x^ - (1 + 'iujf)y^+ (1 + 3a;/) - 3e(l + ujf)z^ + 21VQx\ (29) 

^ = ^[-3(1 - ujs)x^ - (1 + ?>ujf)y^ + 3(1 + ujf) - 3e(l + ujf)z^ - v^7x] (30) 



where 7 = X'y /\/2, C, = \^/ \/2. The various critical points and their properties including the stabilities 
of the above equations are summarized in Table 1 with their eigenvalues being given as follows: 

• point A 

/il = -(w/-l), /i2 = ^, /i3= '- (31) 

• point B 



3(1 - Uf) ^ ^l\-l + a;/)(24 + 24a;/2 + 7^2 + ^^(43 + 9^2)) 
^i'2 - I ± -^2 (32) 

7 + 3a;j7 + 6(1 + Uf)l 
27 

• point C 



3(1 -c^/) , v/e(-l + c^/)(4(l + ^ujff + 3(5 + 27ujf)e) ,.,^ 

/ii,2- ^ ± ^, (34) 

7 + 3u;^7 + 6(l + a;/)^ 

/^3 '■ — 27 ^-^-^^ 
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mi 



/^i 



ointD 


6 + f 


2 


oint E 



/i2 = -3 - 3u;/ - 7 , /i3 



-2 .. 2 + 7^ + 27^ 



(36) 



/ii = -1 - 3u;^ - 6p^ ^2 = l-3f-^, /i3 = -2-3f 



(37) 



point F 



/^l,2 = — ± 



67e(7 + 2O')-37' + 107e + 33e^ 



/W3 



7 + 2^ _ 7 + 2e 

7 + 3w/7 + 6(1 + ujf)S, 



7 + 2e 



(38) 
(39) 



Note that all of the above eigenvalues of the critical points do not have any e-dependence even though 
the various critical points themselves carry its dependence. 

Table 1. The classification and the properties of the critical points. 





X 


y 


Z 


Existence 


Stability 


n^ 


W,j6 


A 











all f and 7 


Saddle point 





- 


Bl 


f (1+^/) 





v^ 


e = -1 


-'T^"<7<0(f>0), 
Stable spiral 


3(l+a,f) 
72 


'^Z 


B2 


l|(l + ^/) 





>K-i) 


e = +l 
UJ[ < -1 


7'>-3(l + t^/), 
Saddle point 


3(l+ojj-) 
7 = 


w/ 


V 272 


C 




V »e 





^/ < "1 


7<-n^fSr (€>0)> stable 


1 + 3^; 


UJJ 


l + 3uf+652 


D 


-^ 





H'^'-^) 


e = +l 


Stable 


1 


-1-f 


E 


«\/l 


^1 + ^ 





all { and 7 


^^ <7 (|'>0), Stable 


1 


-1-2? 


Fl 


2 


/2 + 7(7+25) 
V (7+25)2 


^6(2-35(7+25)) 
3(7+25) 


all 7, e = -l 


-2f < 7 < ^5=|£ (f > 2^2), Stable 


1 


-1 + ^- 


76(7+25) 


F2 


2 


/ 2 + 7(7+25) 
V (7+25)2 


V6(-2 + 35(7+25)) 
3(7 + 25) 


all 7, £ = +1 


f > -v^^ + V;_16+37= (- ^ Q^^ g^^^ig p^^^^ 


1 


1+ "^ . 

^ 37 + 65 


V6(7 + 2e) 



In Table 1, the second to fourth columns denote energy contents of the dark energy sector, the fifth 
column shows the conditions for the existence of the solutions of the critical Equations (28)-(30). The 
critical points B and F are further divided for e = ±1. The sixth column is about the stability condition 



for each eigenvalue. In the seventh column, Vl^ 



r 



ez^ is the density parameter for the 



dark energy, and the last column displays the equation of state for the dark energy. First, we note 
that the results of phases A, B2, D with e = +1 and y = agree with those in the literature [23], 
which correspond to the phantom with a positive potential. In particular, the system approaches the 
scalar-dominated solution D with non-relativistic dark matter [uf = 0). The point F2 also corresponds 
to scalar-dominated solution but with a non-zero value of y. 

Let us focus on the e = — 1 case from here on. Then, points E and Fl both give scalar-dominated 
solutions. Which solution will the system choose depends on the allowed parameters for the existence 
and stability conditions and on the initial conditions. For solution E, z = and cUff, approaches to —1 
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when ^^ — > 1/3. On the other hand, for solution F, all x, y and z contribute to the dark energy and 
co^ converges to —1 when ^ goes to infinity. In the next section, we present an explicit solution of 
Equations (23)-(25), which prolongs from Bl at early times to E at late times. 



5. Exact Cosmological Solution 

In order to investigate the exact solution for Equations (23)-(25), we first assume that a = e'^^'^/y/J, 
which corresponds to H = n^cj). Here / is an arbitrary constant at this stage. In the case of e = —1, 
from the above ansatz one can find an exact solution as follows: 

0(t) = 0(Q) + ^t+ ^ In 

where C is an arbitrary constant and Vq and po are given by 

Vo - 

-ACf 



1 + C'e-3?(i+'^/)t/« 
1 + C 



K^X^{l + ujf) \ -AC 



Po 



3C 



(40) 

(41) 
(42) 



with ^"^ = 2/3A^, 7 = 3(1 + coj)^. It is to be noticed that C should be a negative value to preserve 
the weak energy condition. In this case, introducing a time scale defined by |C| = e~^«^^~'^'^^^^**. 



Equation (40) and a = 

m 

a{t) 



dK?0 



/v7 reduce to 



In ( sinh 



3(1 + a;/) e 



3<(l + w/) 



a sinh 



2 



K 



K 

2 

3(l + LJf) 



fit + u 



m 



(43) 
(44) 



where <^(0) = </)(0) — 2 In sinh( 



'3(l+a;/) g 



v7i*) /(3(1 + a;/X) and ~a = e<~^^^) / ^ 



We note that for t — )■ —t* the solution (44) behaves as a ~ j(2/(3{i+wj))^ which shows that it represents 
an expanding Universe with cosmological fluid ojf, whereas at late times, it represents an accelerating 
Universe with a{t) ~ e^^*/''. Now, let us check that the above solution indeed corresponds to the one 
that prolongs from Bl at early times to E at late times, as was mentioned in the previous section. To this 
end we first substitute the solutions (41)-(44) into the dimensionless quantities x, y, z of Equation (26). 
Then one can find 

1 



X 



v^Ae 

1 sinh[:^i^fv7(t + t)] 



l3{l+uJf) g 



Aecosh[^(^+" 


^!v7(^ + ^= 


.)] 


1 


V l+LOf ' 





v^Aecosh[^(i±^i/7(t + t)] 

When choosing ^ = ■\/2/-\/3A(^ = 1/v^), we see that at early times, i.e., t -^ 
X, y, z of Equations (45)-(47) are given by 



(45) 
(46) 

(47) 
-t*, the quantities 



X 



V2' 



y = 0, z 



1-^/ 
2(1 + ^/) 
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and at late times (t — )■ oo) they become 



v^ ^ V2' 

If we identify 7 = (1 + cij/)-\/3, and substituting ^ = l/VS into Bl and E of Table 1, we exactly find 
the above values of a;, y and z. Therefore, we conclude that our solution corresponds to the phase Bl at 
early times and to the phase E at late times. 

We now determine the allowed range of cuf for this exact solution by comparing with the current 
observational data. It turns out that the non-relativistic dark matter with cuf = Ois, not consistent with the 
current value of to- To carry out this in detail, we recall the density parameter of the scalar field given by 
(we set U ~ 0) 

fi^ = — x^ + y'^ — z'^ 
where B = 3(1 + Uf)^y/f/2K with ^ = \/2/\/3A and the Hubble parameter 

We first comment on the dependence of ^2^, and H of Equations (48) and (49) on the initial conditions. 
The exact solutions. Equations (43) and (44), show the dependence of (p and a on the initial conditions. 
However, the dependence of the quantities x, y, and z on the initial conditions are absent as can be 
seen in the expressions of Equations (45)-(47). This may be due to the fact that our exact cosmological 
solution satisfies the condition a = e^" j \f] and is of a non-perturbative nature. The only remaining 
dependence is through the variable t*, which we have set to zero without loss of generality. Therefore, 
r^^ and R of Equations (48) and (49) are insensitive to the initial conditions as far as the non-perturbative 
solutions (43) and (44) are concerned. This aspect is rather unexpected because there is a priori no reason 
why the quantities x, y and z defined in Equation (26) should not depend on the initial conditions of 
and a. 

In the above two Equations (48) and (49), we treat ^^1^, Hq at t = to as input parameters to 
determine ^^/J and especially cuj. Since the observational data [25] fi^^o — 0.726 ± 0.015, Hq ~ 
2.28 ± 0.04 X 10~^^s~^ and to ^ 4.33 ± 0.04 x lO^^s are given with experimental uncertainty, these 
equations give permitted range of cof. We use the strategy that we first pick up a specific value of cuf. 
Then, Equations (48) and (49) will give allowed ranges of B for fi^ and Hq respectively. We can plot 
these quantities in (B, t)-plane. If there exist intersecting region, then, this value of cUf is allowed. The 
results are displayed in Figure 1. We find that for values of cuf < 0.13 and cuj > 0.22, there does not 
exist a region intersected by the bands of 1^^ 0, Hq and to. Therefore, we conclude that the allowed value 
cof consistent with the current observational data is given by 0.13 < cuf < 0.22, which corresponds to an 
exotic cosmological fluid with non-vanishing pressure. 
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Figure 1. The plot of VtK^Q{hlue) - 0.726 ± 0.015, Ho{green) ~ 2.28 ± 0.04 x lO-^S"^ 
on the (B, t)-plane for uj = 0.13(left), cuf = 0.18(right), uj = 0.22(lower). The red band 
is the current age with uncertainty to — 4.33 ± 0.04 x lO^^s. 




6. Conclusion 



In this paper, we have investigated various cosmological consequences of the de Sitter nonlinear sigma 
model with the exponential potential and cosmological fluid term. It consists of a phantom and triplet 
scalar fields. We have displayed that the presence of the triplet a fields with isotropic spatial dependence 
provides contrasting features to the phantom model without such fields. For example, the energy and 
pressure coming from the spatial variations could prevent the Big Rip singularity from being developed. 
In particular, we found an exact cosmological solution in the case of a negative potential. In this solution, 
the Universe undergoes a power law expansion at early times as in A— CDM. But the difference is that 
unlike the A— CDM we have a non- vanishing scalar energy contribution (x ^ 0, z ^ of Bl phase) 
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even in the matter dominated epoch. At late times, a complete dark energy-dominance is achieved with 
il^ = 1 and u^ = —1 with a suitable choice of the parameter. 

The Hubble parameter Equation (49) is exactly the same as H of A— CDM in Equation (11) when 
we identify the parameter as ^v^37/k = V^. However, with this choice, one can check that the energy 
density pf of cosmological fluid is different from the matter density pm of Equation (11) in A— CDM, 
and is given by pj/pm = 1 + l/('^/ + !)• Note that for < cof < I, Pf is always greater than pm. Also, 
the dark energy density from Equation (48) approaches A/k^ asymptotically from below. These show 
that even if the scale factors in both cases behave exactly the same, there exist qualitative differences 
between the two approaches. 

We also have shown that to be consistent with the observational data, the equation of state for the 
cosmological fluid has to be within the range 0.13 < cuf < 0.22 unlike A— CDM with a dust like matter. 
Thisrangeofu;/- appears in the literature [31,32] corresponding to cosmic strings with —1/3 < coj < 1/3 
or domain walls with —2/3 < cOf < 1/3. We notice that in A— CDM model, the dark matter has the 
perfect fluid form and admits a barotropic equation of state. In this case, the dark matter should be 
pressureless in order to be in accordance with the observational data. However, it is known [33,34] that 
the dark matter can also be described by an anisotropic fluid but only in the case of non-zero effective 
pressure or a polytropic equation of state [35]. On the other hand, our exotic cosmological fluid not only 
satisfies a barotropic equation of state but also takes a perfect fluid form. In spite of this, it should have 
non-zero pressure to be in agreement with the observational data. 

We conclude with the following remark. We found that with cuf = the evolution is exactly that of 
A-CDM. But this does not fit the observational data as was discussed, e.g., cof has to be in the range of 
0.13 and 0.22 to give the observed age of the universe. We speculate that this feature of deviation from 
the A-CDM for general time- varying dark energy density holds in general. That is, in models where dark 
energy density varies, the dark matter with cOm = might have some difficulty in fitting the observations. 
Nevertheless, this does not imply that the dynamical dark energy models with non-zero equation of state 
dark matter must be excluded, and it remains to be seen whether the exotic cosmological fluid considered 
in this work could be related to the realistic dark matter candidates. 
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